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INfRODUllON 

The  design  of  an  open  marine  propultor  11  a  complex 
process.  involving  structural  and  hydrodynamic  conadera- 
( ions  1 1 .  J 1  hoi  the  hydrodynamic  considerations  during  most 
of  the  preliminary  design  process,  approximate  models  of  the 
lifting  surfaces  are  employed, e  g  the  lifting-line  model  (3, 4) 
tor  powering  considerations,  and  two-dimensional  flow  over 
equivalent  Made  sections  for  cavitation  performance  More 
sophisticated  models  of  the  lifting  surfaces  are  used  for  pre» 
dieting  fluctuating  loads!  Si  and  some  cavitation  predic¬ 
tions^!  I  hese  approximate  models  have  been  acceptable 
during  the  preliminary  design  proceu  and  provide  a  bails  for 
choice  of  ihe  maximum  diameter,  advance  coefficient  and 
radial  variations  ol  chord,  skew-angle,  rake,  thickness,  and 
circulation  distribution  I  he  chordwise  variation  m  load  has 
usually  been  selected  during  this  prehminaiv  stage  and  is 
often  based  on  cavitation  and  propulsion  considerations. 


(  hordwisr  velocity  difference 
scaled  to  give  unit  magnitude  when 
integrated  <0  r*»s»  the  c  hord 

\  rmr  bound  In*  rement  to  pitch 
angle  when  radial  inflow  exists 


f-or  the  final  stage  of  the  design,  the  meanline  distribu¬ 
tion  and  radial  pitch  variation  are  determined  corresponding 
to  the  selections  for  load  and  geometry  already  available.  To 
derive  a  geometry  which  accurately  produces  the  specified 
load  distributions,  a  lifting-surface  model  of  the  blades  is 
required 
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Several  procedures  already  exist  for  performing 
lifting-surface  calculations  for  wide-bladed  open  marine  pro- 
pw  Ivors  In  particular,  two  different  appioachet  to  the 
analysis  for  blades  with  arbitrary  locations  in  space  have  been 
presented  by  Ke rwin  (  i )  and  McMahon  ( R)  Kerwin's  numeri¬ 
cal  analysis  procedure  is  based  on  three  fundamental  assump¬ 
tions  1  I )  that  the  continuous  loading  distribution  on  the 
nonplanar  blade  surface  can  be  adequately  approximated  by 
a  multitude  of  discrete  straight  lines  of  constant-vortex 
strength  and  that  the  source  distribution  arising  from  the 
thickness  distribution  can  be  similarly  approximated,  (2)  that 
the  minimum  required  spacing  between  lattice  elements  along 
the  chordlme  is  -  2  degrees,  and  ( 3 1  that  the  resulting 
meanline  shape  for  a  given  chordwise  load  is  similar  to  the 
two-dimensional  shape  for  the  same  chordwise  load  The 
first  two  assumptions  are  not  acceptable  for  very  narrow 
blades  for  »  blade  with  a  20  degree  pitch  -ngle  at  the  0  4 
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radius  and  a  chord  to  diameter  ratio  of  0.0S,  the  2  degree 
spacing  equals  increments  of  about  1/3  chord  length.  The 
last  assumption  permits  calculations  to  be  performed  using 
only  a  few  points  along  the  chord  and  the  two-dimensional 
shape  is  fitted  to  the  data  at  these  points.  The  resulting 
computer  code  is  relatively  quick  running  and  produces  a 
geometry  which,  in  practice,  has  an  overall  speed  and  power¬ 
ing  performance  generally  within  a  few  percent  or  so  of  the 
predicted  values,  with  a  general  tendency  to  produce  a 
greater  thrust  than  predicted.  The  procedure  of  McMahon 
employs  continuous  distributions  for  the  loading  and  thick¬ 
ness  functions  and  calculates  the  meanline  from  the  induced 
velocity  Consequently,  data  at  more  chordwise  points  are 
required  to  define  the  pitch  and  meanline  distributions. 

The  resulting  computer  code  is  lengthy  to  run  but  has  shown 
remarkably  different  meanline  shapes  from  the  two-dimen¬ 
sional  one  at  the  hub  and  tip  region  of  the  blade  where  the 
mean  lines  can  be  s-shaped  (8).  Two  models  were  constructed 
and  experimentally  evaluated  to  provide  data  on  the  relative 
cavitation  and  propulsion  performance  of  designs  having  the 
same  input  specifications  but  final  geometry  according  to 
the  Kerwin  and  McMahon  procedures.  Some  inconsistencies 
occurred  in  the  experimental  measurements  but  the  thnist 
was  closer  to  the  predicted  value  and  the  operating  point 
centered  in  the  cavitation  bucket  for  the  model  designed  by 
the  McMahon  method  Hence,  the  determination  of  specific 
meanline  and  pitch  distributions,  instead  of  fitting  the  two- 
dimensional  meanline,  is  considered  to  he  a  superior  pro¬ 
cedure  when  the  design  is  based  on  a  narrow  range  of  permis¬ 
sible  operating  conditions  and  the  delay  of  cavitation  is 
critical. 

Because  the  numencaJ-anaJysis  procedure  employed 
by  McMahon  results  in  lengthy  computer  runs  and  Kerwin’s 
procedure  is  not  acceptable  for  narrow  blades,  alternative 
numerical-  ana  lysis  schemes  are  investigated  in  this  paper. 

In  addition,  a  detailed  description  of  the  flow  field  across 
the  blade  surface  was  desired  as  input  into  boundary-layer 
calculations.  Two  different  numerical  analysis  schemes  are 
described,  each  involving  an  expansion  of  the  singular  kernel 
about  the  singular  point  Both  approaches  employ  integra¬ 
tion  of  the  specified  thickness  slope  and  load  distribution 
over  the  reference  blade  in  the  radial  direction  first  and  the 
remaining  chordwise  integration  then  takes  the  form  of  the 
velocity  component  corresponding  to  two-dimensional  flow 
modified  by  the  presence  of  an  induction  factor  in  the 
integral  Regular  integration  techniques  are  employed  for 
the  other  blades  and  the  shed  vortex  sheet  The  induced 
velocity  components  are  appropriately  combir-d  and  inte¬ 
grated  to  obtain  meanline  shapes 

The  present  investigation  describes  the  real- fluid 
flow  about  a  rotating  system  of  lifting  surf  aces  having  both 
loading  and  thickness  Several  approximations  are  made 
The  first  of  these  is  the  mathematical  model  for  which 
potential  flow  equations  arc  employed  and  the  solution  to 
first-order  in  t hie  knesvto-ch  >rd  ratio,  camber-t»>-chord  ratio 
and  difference  in  pitch  and  flow  angles  derived  (  ompari 
sons  with  experimental  results  for  other  lifting-surface 
configurations  lead  to  confidence  in  this  linearised  approx¬ 
imation  In  addition  to  this  mathematical  model,  further 
approximations  occur  m  the  numerical  analysts  C  onfidence 
in  the  numerical  analysis  procedures  is  luslified  by  cornpar 
ison  with  analytical  solutions  or  experimental  results  1  hat 
is.  results  are  sought  from  some  disc  retired  numerical- 
analysis  procedures  involving  N  by  M  approxima  ions,  which 
have  converged  to  within  some  specified  tolerance,  t .  of  the 
real  or  analytical  value  of  fhe  quantity  investigated  Math 
ematuaily  this  may  be  stated 

'  f  t  x.  y  1  -  f^  M  t  x  v  i  *  t 


{(x,  y)  on  the  surface  S 
N  >  N(| 

M>MC 

where  fN  M  =  the  approximate  calculation  of  a  particular 
quantity  f 

S  =  a  region  of  the  surface  of  interest 

N0,  M()  =  minimum  numbers  of  the  discrete  approx¬ 
imations  for  which  the  computed  results  are 
within  €  of  the  values  for  f 

For  rotating  lifting  surfaces,  neither  measured  nor  analytical 
solutions  exist  for  details  of  the  flow  field  on  the  blade 
Hence,  comparisons  will  be  made  with  other  procedures 
It  is  assumed  that  numerical  solutions  which  employ  in¬ 
creasingly  greater  pointwise  definition  of  the  input  variable 
without  change  in  computed  values  have  converged  and  that 
the  solution  has  converged  when  a  smooth  curve  can  he  drawn 
through  point  values  in  both  the  cnordwise  and  radial  direc¬ 
tions.  These  assumptions  are  believed  to  be  necessary  but 
not  sufficient  for  convergence. 

In  the  following  sections,  the  mathematical  model 
of  the  flow  field  on  the  blade  surface  is  first  reviewed  and 
numerical-analysis  techniques  for  evaluating  both  regular  and 
singular  integrals  are  described.  A  FORTRAN  computer  code 
is  discussed  and  sample  calculations  using  this  code  are  pre¬ 
sented.  From  example  calculations,  it  is  found  that  greater 
accuracy  in  the  integral  evaluations  is  required  for  the  deter¬ 
mination  of  smooth  pressure  distribution  curves  than  for  the 
shape  of  the  meanline  and  the  pitch  distributions.  The  choice 
of  a  particular  chordwise  loading  distribution  is  shown  to 
have  an  effect  on  the  meanline  shape  and  the  pressure  dis¬ 
tribution.  The  effects  of  rake  and  skew  are  shown  to  be 
important  on  both  pressure  distribution  and  meanline  shape 
A  particular  thickness  function  has  hardly  any  effect  on  pitch 
or  meanline  but  a  significant  effect  on  pressure  distribution 

MATHEMATICAL  MODEL  THICK  LIFTING  BLADE 

The  mathematical  model  of  a  system  of  rotating  lifting 
surfaces  advancing  in  an  unbounded  irrotational  flow  field 
with  an  inviscid  fluid  has  been  developed  on  a  formal  mathe¬ 
matical  basis  by  Brockett  19)  A  reformulation  of  that  analysis 
in  terms  of  non-dimensional  surface  coordinates  is  presented 
herein  for  completeness  The  propulsor  is  assumed  to  be  ade¬ 
quately  represented  by  the  blades  alone,  i  e  .  neither  fhe  hub 
nor  fillet  from  the  blades  to  the  hub  is  included  in  the  blade 
specification  The  onset  flow  is  assumed  to  he  directed  along 
the  axis  of  rotation  hut  a  new  feature  included  herein  is  that 
it  may  have  a  small  radial  component  Overall  geometry  nota¬ 
tion  generally  follows  the  definitions  given  in  Reference  10 

Coordinate  systems  are  constructed  with  the  same 
orientation  as  in  Reference  and  in  particular,  the  helical 
coordinate  system  <fc| ,  r)  rotating  with  the  blades  is 
shown  in  Eigure  1  Unit  base  vectors  in  a  right  handed 
(  artesian  reference  frame  are  the  customary  t  i,  j.  kl  where  i 
is  along  the  x  axis  and  is  positive  pointing  aft.  i  is  along  the 
v  axis  and  k  is  along  the  /  axis  which  is  generally  along  the 
reference  blade  Unit  vectors  along  the  helical  coordinates  are 

e ,  r  sin  i  ♦  cns<0p  i  i  » 

v  ^  -  c  <»s  0p  i  ♦  sin  0p  '  i 

r  f  -  sin  i  ♦  *  os  it  k  t  t » 

where 
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VERTICAL  BLADE  REFERENCE 

*  /  A  LINE 

*  AFT 


tl  »  c|xc-0.5) 


l-ig  1  Lifting-surface  geometry 

The  blade  surface  is  given  by 

=  F («, . r>  t: 

=  Ec  (£, ,  rl  t  hT  (£, ,  r)  (f 

where 

Fc  is  the  mcanline  shape,  and 
Ft  is  the  thickness  shape 

In  the  analysis,  it  is  convenient  to  change  the  variables  of 
integration  to  (  x£.  xR  1  instead  of  ({, .  rl.  where 

1 1  =  c  ( x  -  0. 5 »  V 


r  =  Dx„/2 

c  =  chordlength  at  radius  r 
D  =  maximum  rotor  diameter 


The  position  vector  of  a  point  on  the  hlade  surface 
described  by  Equation  f  5 )  is 

'  =  D{[d  *  D  <\ -05lstn*p-  bp  cosOpji 


and  a  normal,  directed  out  from  the  blade  surface,  is  (<#.  II) 
N  *  t  x  • 

d  XL  d  X|J 

where  the  plus  sign  is  used  tor  the  suction  side  of  the  blade 
and  the  negati  e  sign  for  the  pressure  side  of  the  blade 
After  some  effort  it  tan  be  shown  that 


c  a  i  /I)  d  c/D  a  f  /D 

Nn  -  -  2 - +  2(x  *0.5)  - - - - 

K  I)  d  x^  c  d  Xfl  3  xt 

llil/l)  Tc  a  l-:/l> 

♦  2  - —  -  cos  ip p  +  - -  sin 

dxR  Ln  8\  J 

T/cv2  id  I./D1 

[Vl)/  1  l)ftxL 

[d  P/I)  cos2  ftp  sin  pp  cos  </>p 
dxR  »xR  xR 

dfls  i  < xt  -0,5 )  cos  <Pp  sin  Up 

-  Xt)  “  2  " 

d  xR  xR 

c  a  t/D 

•  —  sin  -  - -  cos  ®p 

D  r  n  xt  r 
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constant-pitch  blade  which  is  neither  raked  nor  skewed 
In  Pquations  ( 1 0)  and  (11) 

iT  =  the  total  rake 

P  =  the  pitch  of  the  hlade 

0p  =  the  pitch  angle.  0p  =  tan*1  I P'(wxRDu 

8  =  the  angular  position  of  a  point  on  the  hlade 
surface,  a  function  of  both  x  and  xR 
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7  *  [r>  T 

*77',n  «p  /r 

=  the  skew  angle,  a  function  of  Xp 
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where  the  symbol ^  means  symmetry  resin  1 1 ■  »»* 
the  limiting  region  whkh  excludes  the  singular  it  \ 
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This  known  value  a(  the  singular  point  allows  a  straightfor¬ 
ward  analysis  procedure  to  be  undertaken  using  the  pro*. cduies 
previously  described 

Some  convergence  problems  near  the  leading  and 
trailing  edges  and  over  much  of  the  surface  for  narrow  blades 
(maximum  c/D  ^  0  05)  have  been  resolved  by  computing 
the  hnean/ed  form  of  F  (Equation  100)  over  the  entire  blade 
and  adding  a  correction  term  which  is  the  difference  between 
the  actual  integrand  and  this  linear  approximation  This 
option  has  been  included  in  the  computer  program  and  is 
defined  as  “linear  approximation-plus-difference  “  When 
conventional  integration  techniques  are  used  everywhere 
except  at  the  singular  point,  where  Fq  nation  ( 103)  is  re¬ 
quired.  the  procedure  is  defined  as  “direct  " 

For  the  trailing-vortex  sheet,  a  regular  integration  can 
be  performed  since  no  singular  points  occur  on  the  sheet 
The  strength  of  the  vorticity  is  given  by  Equation  (  5^1  and 
the  induced  velocity  field  is  given  by 
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Sp.  »  .1  <  "i»i|  'it.  1  1  1 . r « •  t'i  t-.iitf*  •  in  t  ,  |,i  r  i  p|  .  1 

1  rsult  nif  -ii  1  tti.i « imniti  ti.if »?,  if  1/  ■  .-it  ^  4i  \  II  i  f  .  -  :  ,<■ 
I'ffs,  nti  -  I  I  f.  *.)n.  1  i*  1  .it  op,  .ill  ;  -t  1  1  ,t  t  -  l  \  "  s  It-  aiM 
ii>.  nit  "fils  1  ifti  p*  1  fnl  -)|t  1  •  1 1  M,  r  I-  pit  I  I  -*ri  t  *  f  *  .il.i 

•*h"ut  I  f  ■  1-  s.iiiii  ■!'  t  -  mi  ii.  I  I  <  .1  k'  I  uv  1 1  -  n  ,  itii  1,.  il  n  1 1 »  si  v 

I  hi  •  -  input  1  . 1  pit  fi  fi-mn,  1  1 .  1  t  vfc  p,  t  .  t'l  U-  vs  r  t  ,|i 
.  "Miplltril  t>\  kl  tXUM  v  MHth..,|  Sin- (  illl-i.t  h\hf.l  i- np.*r, 

•  ik  *  .ii  I )  1  VSN  I H  ti.  .  1 .1 1 ,  has  h*  ..-fi  that  k  1  *  m  pi  - .  ■  1  ,r  < 
pf  kIn.  t-x  Jpsifpis  t fiat  at*  »frn«  r.il|\  s|t^hf|v  .  ,*«  rpit  in  J 

pt  fh.tps  sum*-  iffjpmv, m,  nf  m  prf f <*rrn.ifi  .  ma\  t*r  .  vp«  1  ! 

living  the  pn  s«-nt  iiifth"*! 


PrrJl,  tpifis  .if  till-  n»f  h  -M>'l  .ill  1 1  I  1  *  v  t*u  1«v- 
pf".  i-iluffs.  .U‘v*  l"pt  J  Nti  •  •  ■  in p u 1 1  n p  th*  im-I.i.  c«l  v.  I-  it  -, 
f  i«a  Ut  »m  thr  hi. 1J1  suita.i  yv  in.  h  .  .  *n  t  asm  tlt,  tp'IJ  p<  <ipt 
Jiir.t  .nut  .ippfi'Ximatr  pin' Jif  ti  fi'fp  ,  <  r  <  sli-wi  ’ 

I  jNi  I  l"  i-<  math  fin  s.inn  H"t*i\i  t  it  *■..*'  t-r-  n  t  pi  : 
that  "Vi  tjll  tlii  apprnxmiati  plus  JiMimpm.  1  pt".  p-  v 

ptclctiihlr  vkhvn  Ji  tisc  .  Imrilvusi*  spj,  m  fi"M-n  * . 

S\  I ’/ *  i  u  nartnviw  t'laJps « fii.txittmm  It  iM^'.it. 
inviilvri)  In  this**  situatiuns  thp  Jir*Kl  pt"  i  ilnti’ pr  »Jn . «•' 
lnv.ill\  iTt.itii  values. J  the  uukueJ  vrl<.,  1 1 \  hiKausi  nt  the 
Je.  tease. t  v|'.i.mj!  between  aJi.Kent  lines  of  mtei»rjUi»n  witii  j 
.  nrresp<ifijifi(?  laiK  "I  Juuuiv  in  the  nuni.  ru  al  intent jt  1. nu 
tut  the  resulting  tieat  .insular  integrals  I  his  rttevt  is  illus 
trateJ  in  I  ifiure  ^  whu  h  sh*<ws  values  "t  .me  "I  (lie  helkai 
vomp.ments  .»(  r/u-  average  imfuv<sJ  wh»,  itv  ./ f  tile  '/  t4t- 
taJius  .J  I  lit-  feteteike  hlaJe  lh|svel".itv  ,  *»mp*>nent  is  .lu 
tu  **nk  U i.ul  1  n>*  on  the  blade  itsell  the  etUkls  tiu.kness 
the  other  blades  and  the  shed  vortex  sheet  m  1  MklmleJ 
XII  data  shown  in  subsequent  Inures  ha'«  t'*  en  .oinputed  *'v 
the  'appioxinute  pluvditterenve "  pf»Kedure  although  onlv 
the  ptessure  disltibution  neat  the  leading  ed^e  in  the  tip 
teirinn  »d  the  tdade  was  signif  K  ant  lv  dtlterent  between  Un- 
two  pFondiiies 

( iver.ill  run  time  varies  with  number  »>l  points  number 
of  i'l.tiles  and  blade  width  Silue  omiputei  iisa^e  1  hat  lies  ,»fp 
so  low  the  I  x  I  x  I '/ atfav  M/e  is  ti\  oinuiended  I  or  a  n.iriow 
blade  I  lie  linear  '  approximation  plus  difference  procedure  •' 
teMMiuneriiled  ami  file  run  lime  mav  iiu  tease  b\  a  lew  turn 
dred  seinnds  b\\ jtixe  ol  speviaKare  taken  with  the  shut 
vortex  sheet  .ahulations  (  oinpwterexe.  ition  tune  lot 
kerwm  s  program  iv  unknown  for  rhe  Hur/ou^hs  'DO  hivh 
speid  vompuier  but  |\  estimated  t«»  in-  about  l  Mi  seconds 
1**r  data  .  aU  illations  at  4  ,  hordwise  points  at  H  tadial  stations 
l  of  the  results  shown  in  I  able  I  data  ate  vOtnpufeit  at  I  1 
tadial  stations  with  either  *  I  I  01  I  "  .  hold  wise  points 
ilrpen  Imp'  on  input  data  spe.  itn  atmn 


I  ufther  details  of  the  ^eornettv  *»l  this  1  xample  an- 
ffiveu  m  I  able  II  Radial  v  atiahles  ate  titled  .u.oidme  t‘>  the 
s\mb.»|s  sulip'ested  111  RetereUve  10 
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IHM  I  XSIOS  (>t  \  \ AMP!  !  (  OMPl  I  A  HONS 

In  f  His  section  tt.o  c onsequences  nt  Ihr 

•lrM|£ni  r  nughl  make  both  lor  overall  geometry  am]  lor  Ihr 
hordwise  vjrulion  <>l  Ihr  thickness  distribution  and  loading 
distribution  are  examined  Some  ininmnn  variations  in  Ihr 
location  of  ihr  blade  mid  chord  line  are  investigated  In  defer 
mint'  thr  effects  ■>!  ovrrall  geometry  on  pilch  camber 
pressure  distribution  and  second -order  performance  nrfli 
u  nis  I  hr  variations  arr  unskewed  skewed  and  warped  I 'h 
Hades  wlh  other  input  specifications  the  same  Skewed 
Hades  have  blade  sections  displaced  along  the  pitch  helm 
and  warped  Hades  fuse  Made  sections  displaced  ctrcum- 
terent tally  in  Ihr  plane  at  x  '  0 


Inpul  quantities  and  selected  nutpul  are  shnwn  in 
I  able  II  tm  a  warped  blade  similar  In  NSR|M  Model  44^H 
cam!  one  similar  In  ihe  example  nl  Reference  ’)  I  <>r  an 
unskewed  blade  I  tie  mluinn  labeled  I  I  I  S  the  skew  angle 
Wv  would  hr  /ern  anti  I  nr  a  skewed  blade  Ihe  column 
labeled  KAK  I  l>  the  total  rake  ij  I)  would  be  equal  to 
P  •  I  '#l)i  In  I  igure  4  the  computed  pilch  and 
,  amber  rjlms  are  shown  tor  these  various  overall  geometries 
with  all  other  input  Ihe  same  js  in  I  able  II  In  figure  S 

m.  i!  is  ..|  itii  vh<>i<lwisr  load  disliihiitii -n  ltiul  hordwis. 

thk  kness  furu  lion  on  pdc  h  and  c  amber  are  shown  I  he 
ellevt  nt  rake  ami  skew  on  pilch  and  «.  amher  follows  known 
trends  <  /4»  I  he  efte»  t  ol  I hK kness  disTnhution  on  pile  h 

and  camber  is  negligible  and  ihe  effect  ol  elliptic  loading  is 
I"  reduce  Ihe  pitch  and  increase  the  camber  a.  would  he  Ihe 
vase  m  two  dimensional  flow  at  the  ideal  angle  lot  a  given 
lilt  coefficient  In  I  igure  h  Ihe  pile  h  and  camber  change  is 
shown  tor  another  modilicalion  ol  Ihe  warped  blade  Since 
.i  large  change  occurs  in  the  pitch  from  rhr  input  spec  ific  a 
lion  t  f  able  fl  i  In  Ihe  computed  values  ( f  iguo*  t»  cornputa 
lions  were  performed  with  Ihe  singularities  distributed  on 
the  Made  reference  surface  at  a  pilch  taken  from  figure  4 
I  his  change  in  pitch  place-,  ihe  singularities  nearer  Ihe  final 
blade  surface  lo  have  unitormilv  in  the  calculations  the 
pilch  angk*  ol  Ihe  shed  vortex  sheet  was  taken  as  d  the 
advance  angk*  ol  the  shed  vortex  sheet  iln  t  igure  4  ihe 
shed  vortex  sheet  was  taken  lo  be  al  Ihe  input  pilch,  which 
is  d,  Ihe  pitch  angle  derived  from  the  solution  ol  a  straight 
radial  Idling  line  representing  each  blade  »  1  he  ,  hange  in 
pitch  angle  of  Ihe  shed  vortex  wake  from  d,  to  d  produces 
a  slight  iru  tease  ol  pilch  near  Ihe  hub  I  compare  data  in 
I  igures  4  and  <>»  A  change  in  the  pilch  »»l  Ihe  blade  reference 
surljie  to  Ihe  values  shown  by  the  dashed  curve  hi  figure  4 
produces  a  significant  reduction  in  computed  pile  hand  a 
compensating  increase  in  camber  near  the  root  with  negh 
giblv  change  in  either  pitch  » »r  camber  trom  about  x#  0  * 
to  the  tip  Hence  the  orientation  of  the  tree  vortex  sheet 
and  blade  reference  surface  have  significant  ette,  Is  on  Ihr 
pilch  and  camber  values  onlv  near  I  tie  hub 
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It  is  often  convenient  to  have  an  orthogonal  coordi 
nate  system  on  the  surlace  ol  the  blade  In  particular  tor 
performing  boundary  layer  computations  an  orthogonal 
coordinate  system  with  one  variable  along  the  streamlines 
reduces  the  number  of  terms  in  t he  governing  equations 
fo  determine  the  differencial  equation  of  the  streamline 
path,  lei 

\K  =  vl\>  M 

be  the  radius  of  the  streamlines  as  a  (unction  ol  the  c hord 
wise  coordinate  Then 

s*i\j=sixc^(\h  (I'li 

is  the  position  vector  ol  the  streamlines  on  the  Made  surface 
Hence  a  tangent  to  the  streamline  is 
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I  t  this  tangent  vector  to  he  parallel  to  the  velocity  vector 
on  the  surlace.  the  vector  cross  product.  M  X  q.  must 
be  /cru  Hence.  Jor  the  velocity  on  the  blade  surface  given 
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Thus  the  slope  ol  lines  on  the  surtav 
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One  now  has  differential  equations  to  determine  .in  oitlio 
onal  network  over  the  Made  surlj.e  I  he  dillerenli.il  .11 . 
length  along  the  streamlines  is 
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